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Abstract. Precision measurements of nucleons provide constraints on the Standard Model and can discern the signatures 
predicted for particles beyond the Standard Model (BSM). Knowing the Standard Model inputs to nucleon matrix elements 
will be necessary to constrain the couplings of dark matter candidates such as the neutralino, to relate the neutron electric 
dipole moment to the CP-violating theta parameter, or to search for new TeV-scale particles though non-V— A interactions in 
neutron beta decay. However, these matrix elements derive from the properties of quantum chromodynamics (QCD) at low 
energies, where perturbative treatments fail. Using lattice gauge theory, we can nonperturbatively calculate the QCD path 
integral on a supercomputer. In this proceeding, I will discuss a few representative areas in which lattice QCD (LQCD) can 
contribute to the search for BSM physics, emphasizing suppressed operators in neutron decay, and outline prospects for future 
development. 
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The Standard Model is one of the great successes of modern physics, and all experiments done so far are consistent 
with the predictions of the Standard Model (SM) within reasonable errors. However, there are certain deficiencies of 
the Standard Model that lead us to consider physics beyond the Standard Model (BSM). The LHC has been built to 
discover the mechanism for electroweak symmetry breaking and to probe certain potential new BSM particles which 
appear at the TeV scale, such as some candidates for dark matter. Such high-energy colliders in general are well suited 
to production of high-mass particles directly at their resonances and for probing processes with cross-sections that 
scale with energy. However, there are many low-energy precision experiments that also have the ability to probe BSM 
physics. These generally fall into two classes: experiments that can very precisely measure physics that is precisely 
predicted by the SM, such as muon g — 2, the proton weak charge, and flavor-changing processes involving the CKM 
matrix. Another class of experiments can look for tiny signals in places that the SM says should be either undetectably 
small or exactly zero, such as the neutron electric dipole moment, or non-V— A contributions to neutron beta decay. 
These low-energy experiments, which usually involves detections through interactions with nuclear matter, will involve 
the low-energy nonperturbative properties of QCD. Lattice QCD is an ideal theoretical tool to calculate these from the 
first principles of the Standard Model. In this proceeding, we give three examples where LQCD can contribute to the 
search for BSM physics using nucleons. 

The Strange Proton: Recent evidence suggests that the composition of dark matter must be cold, with velocities 
that are non-relativistic, giving it a Compton wavelength above the width of a proton. One promising candidate for 
dark matter is the supersymmetric particle called the neutralino. If dark matter is neutralinos, they will interact with 
nuclei at low-energy. A typical value for the cross-section in the Constrained Minimal Supersymmetry Standard Model 
is about 1 zeptobarn (1CT 49 m 2 ). The strange-quark contribution to the proton scalar density and spin are important 
inputs for spin-independent and -dependent cross-sections. In this proceeding, we will focus on the LQCD calculation 
of a s = m s (N\ss\N). 

One can directly calculate the matrix element which involves strange-quark loops and proton correlators. Unfor- 
tunately, such a direct calculation is noisy, so various techniques have been developed to improve the signal (e.g. 
JLQCD, QCDSF, BU, Engelhardt)[l]. Alternatively, one could use the Feynman-Hellman Theorem by studying how 
the nucleon mass varies with the strange-quark mass, taking the numerical derivative dM^j / dm s either by direct SU(3) 
fitting to the baryon masses (e.g. Young/Thomas) or by reweighting strange part of action (e.g. Jung, MILC)[2]. The 
values of a s from various calculations are shown on the left-hand side of Fig. 1 . We perform a global fit to current 
dynamical lattice data (weighted by lattice spacing, lightest m K , dynamical strange and other quality factors) and ob- 
tain (7 S LQCD = 47(8) MeV, which is consistent with the value obtained using only two lattice inputs in Ref. [3]. The 
cross-section constraints from various models in the CMSSM can be found in Ref. [3]. 
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FIGURE 1. (left) Summary of dynamical LQCD calculations of the nucleon sigma term m s {ss) in MeV. The (pink) band indicates 
the global average after weighting individual points according to their systematics. (right) Summary of the latest dynamical d„ 
calculations as a function of m^. The band takes all the lattice points with weighted errobar and globally extrapolates the various 
calculations to the physical pion mass point (the leftmost star point). 
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Neutron Electric Dipole Moment: The neutron electric dipole moment is a measure of the distribution of positive 
and negative charge inside the neutron. To generate a finite nEDM, one needs processes that violate CP-symmetry, 
for example by adding a CP-odd 0-term to the lagrangian. The SM value of this quantity is very small, 10~ 30 e • cm. 
Although experiments do not have the necessary precision to measure the SM value; many BSM models predict 
values that are higher than the experimental upper bounds, including some parts of the parameter space for certain 
SUSY models. 

There are three main approaches to extracting the nEDM using LQCD. One is a direct measurement, by studying 
the electromagnetic form factor F3 under the new QCD lagrangian including the CP-odd term (RBC, J/E, CP-PACS 
2005). This requires us to extrapolate the form factors to q 2 = 0, which can introduce systematic error and exacerbate 
the statistical error. Another method is through indirect measurement by looking at the energy difference for two 
different spin states of the nucleon at zero momentum under an external static and uniform electric field (CP-PACS 
2008). Or finally, one can measure the product of the anomalous magnetic moment of neutron F2(q 2 = 0) and tan(2a) 
(QCDSF) in a CP-violating system. The collected dynamical LQCD results [4] on d„ are shown on the right-hand side 
of Fig. 1; when combining all data and extrapolating to the physical pion mass, we obtain d^ CD = —0.015(5) Qe ■ fm. 
More updated and precise calculations from various groups are currently in progress. 

Non-V— A Interactions in Neutron Beta Decay: A third example of probing BSM with nucleons is to look for 
signatures of contributions due to scalar or tensor interactions in neutron beta decay[5]. The notion is quite similar 
to how Fermi theory led to the discovery of the electroweak interaction and its bosons. Beta decay was originally 
explained by Fermi by adding a new term to the fundamental lagrangian describing 4-fermion interactions. Such a 
theory introduces a coupling, the Fermi constant, that has units of inverse-energy squared, with an energy scale around 
100 GeV. As it turns out, the Fermi theory is a low-energy effective theory approximating the electroweak theory, 
which has 3 vector bosons, the W's and Z. The theory was later probed in high-energy proton-antiproton experiments 
at CERN, and the new particles were found with resonances around the scale predicted by this interpretation of the 
Fermi theory. We can imagine that new particles beyond the Standard Model can be predicted in just such a way: by 
first detecting the low-energy effective operators and later directly probing them in a high-energy experiment. 

The new BSM physics will enter the neutron beta-decay hamiltonian, e? SM 0. ept x of 3 * , with the coefficients 
e = v 2 /(Af SM ) 2 , where Af SM are related to the TeV scale of the particles. The leptonic part is understandable using 
analytic techniques, but the quark operator in the context of the nucleon will introduce some unknown coupling con- 
stants, gj — {p\uOn V d\ri) and g$ = (p\ud\n), which are nonperturbative functions of the nucleon structure, described 
in the SM by QCD. Any deviation from the SM V— A current coming from the new scalar and tensor interactions in 
the effective theory will require knowledge of the couplings g$ and gj to understand. 

The tensor charge (#7-) is also known as the zeroth moment of transversity; experimentally, one can extract the 
contribution from individual quarks as a function of the quark momentum fraction x, at certain Q 2 through processes 
such as semi-inclusive deep inelastic scattering. To obtain gr, one need to integrate over all possible x, that is, from 
to 1, through an interpolation of the few available x at each Q 2 . After combining a few experiments, the best 
experimental estimate (using model inputs) for this quantity is obtained at Q 2 = 0.8 GeV 2 (instead of Q 2 = as in 
the correct definition of gj): 0.77+2'^ with combined uncertainty around 25%. There are also theoretical models that 
attempt to make estimates of the tensor charge, such as Nambu-Jona-Lasinio model and the chiral quark soliton model; 




FIGURE 2. (left) Global analysis of all Nf = 2 + 1 lattice calculations of gj and gj. The left-most points are the extrapolated 
values at the physical pion mass. The two bands show different extrapolations in the pion mass: with and without the lighter PNDME 
collaboration points; the final uncertainty is significantly improved due to the chiral extrapolation, (right) £$-£7- allowed parameter 
region using different experimental and theoretical inputs. 



unfortunately, they are not consistent with each other. And using QCD sum rules, still yields a large uncertainty. The 
other important ingredient for the new interactions is the scalar charge g$. There are no experimental measurements 
of this quantity, and the theoretical estimations (from different model approximations) give rather loose bounds to the 
quantity: 0.25 < gs < 1. 

Lattice calculation of gj and g$ is rather straightforward. The tensor charge has been studied a few times in the 
past using Nf = 2+1 dynamical ensembles in LQCD, while gs has not been studied until now; a summarized LQCD 
results are shown in the left-hand side of Fig. 2. All of these results, by different collaborations, using different fermions 
actions, are generally in good agreement, although the error bars shown here contain only statistical error. We further 
globally analyze all the available lattice data, including the lightest available from PNDME, around 220 MeV, and 
obtain g^ CD — 0.95(5) and g!^ CD = 0.69(9) through chiral extrapolation to the physical pion mass. 

Finally, we can combine the low-energy tensor and scalar charges with experimental data to determine the allowed 
region of the epsilon parameter space. We can combine the current knowledge of g$j and the existing experimental 
data for the + — > + transition and nuclear beta decay (such as j3 symmetry in Gamow-Teller 60 Co, longitudinal 
polarization ratio between Fermi and Gamow-Teller transitions in 114 In, positron polarization in polarized 107 In and 
beta-neutrino ccorrelation parameters in nuclear transitions). Using the g$.r from the model estimations and combining 
with the existing nuclear experimental data, we get the constraints shown the outermost band on the lower-right of 
Fig. 2. There is an on-going ultra-cold neutron experiment studying neutron beta decay at LANL to look for deviations 
from the SM in the Fierz term and the neutrino asymmetry parameter to the level of 10~ 3 by 2013. Combining those 
expected data and existing measurements, and again, using the model inputs of gsj, we see the uncertainties in £57- 
are significantly improve. This shows the importance of the precision experimental inputs. Finally, using our present 
LQCD values of the scalar and tensor charges, combined with the expected 2013 precision of experimental bounds we 
found the constraints on Esj are further greatly improved, shown as the innermost region. This corresponds to lower 
bounds for the scale A$j at 2.9 and 1 1.2 TeV, respectively, for new physics in these channels. 
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